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ABSTRACT 

We report on an accounting of the star formation and accretion driven energetics of 24 fim detected 
sources in the Great Observatories Origins Deep Survey (GOODS) North field. For sources having 
infrared (IR; 8 — 1000 /im) luminosities > 3 x 10^^ Lq when derived by fitting local spectral energy 
distributions (SEDs) to 24 fim photometry alone, we find these IR luminosity estimates to be a factor 
of —4 times larger than those estimated when the SED fitting includes additional 16 and 70 fim data 
(and in some cases mid-infrared spectroscopy and 850 /im data). This discrepancy arises from the 
fact that high luminosity sources at z 3> appear to have far- to mid-infrared ratios, as well as 
aromatic feature equivalent widths, typical of lower luminosity galaxies in the local Universe. Using 
our improved estimates for IR luminosity and AGN contributions, we investigate the evolution of 
the IR luminosity density versus redshift arising from star formation and AGN processes alone. We 
find that, within the uncertainties, the total star formation driven IR luminosity density is constant 
between 1.15 ^ z < 2.35, although our results suggest a slightly larger value at z > 2. AGN appear 
to account for ^18% of the total IR luminosity density integrated between < z < 2.35, contributing 
<25% at each epoch. Luminous infrared galaxies (LIRGs; 10^^ < Lir < 10^^ Lq) appear to dominate 
the star formation rate (SFR) density along with normal star-forming galaxies (Lir < 10^^ Lq) 
between 0.6 ^ z < 1.15. Once beyond z > 2, the contribution from ultraluminous infrared galaxies 
(ULIRGs; Lir > IQ-'^^ Lq) becomes comparable with that of LIRGs. Using our improved IR luminosity 
estimates, we find existing calibrations for UV extinction corrections based on measurements of the 
UV spectral slope typically overcorrect UV luminosities by a factor of ~ 2, on average, for our sample of 
24 ^m-selected sources; accordingly we have derived a new UV extinction correction more appropriate 
for our sample. 

Subject headings: galaxies: evolution - infrared: galaxies - ultraviolet: galaxies - radio continuum: 
galaxies 



1. INTRODUCTION 

Precisely quantifying the stellar mass assembly over 
cosmic time is critical for understanding galaxy forma- 
tion and evolution. While the first studies which set 
out to characterize the cosmic star formation history had 
to rely on observations of the emergent (rest-frame) UV 
emission, requiring significant corrections for extinction 
by dust, the advent of space-based infrared observato- 
ries have allowed us to measure the dust-obscured star 
formation activity directly. The resolution of the cosmic 
infrared background using ISO (e.g., Elbaz et al. 2002), 
Spitzer (e.g.. Dole et al. 2006), and now Herschel (e.g., 
Berta et al. 2010), has uncovered populations of dusty, 
starbursting galaxies with infrared (IR; 8 — 1000) lumi- 
nosities ranging between IQ^'^ Lq < Lir < lO^'^ Lq. 
These cosmologically important luminous (lO^^ Lq < 
Lir < 10^^ Lq) and ultraluminous (Lir > 10^^ Lq) in- 
frared galaxies, LIRGs and ULIRGs respectively, appear 
to dominate the star formation rate (SFR) density in the 
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Universe between redshifts of 1 < z < 3 (e.g., Chary & 
Elbaz 2001; Le Floc'h et al. 2005; Caputi et al. 2007). 

Many of these studies rely on deep 24 fim imaging and 
the use of local spectral energy distribution (SED) li- 
braries to extrapolate the evolution of the total infrared 
luminosity with redshift. In doing so, large (factors of 
~10) bolometric corrections are invoked to scale up the 
rest-frame mid-infrared flux densities and estimate the 
total SFR per unit co-moving volume. Using the rest- 
frame mid-infrared wavelengths, speciflcally near ^^8 /im, 
can be problematic due to the presence of broad spec- 
tral features attributed to polycyclic aromatic hydrocar- 
bon (PAH) molecules. While the rest-frame 8 /xm emis- 
sion correlates with the total IR luminosity of galaxies in 
the local Universe, albeit with a large amount of scatter 
(e.g.. Dale et al. 2005; Smith et al. 2007; Armus et 
al. 2007) and systematic departures for low-metallicity 
systems (e.g., Engelbracht et al. 2005; Madden et al. 
2006), there are hints that this correlation may break- 
down at increasing redshift. In a recent study. Murphy 
et al. (2009) has shown that galaxies in the redshift 
range between 1.4 ^ z < 2.6 exhibit PAH equivalent 
widths that are large compared to local galaxies of sim- 
ilar luminosity. This translates into an overestimate of 
IR luminosities from rest-frame mid-infrared photometry 
which has been noted for subsets of z — 2 galaxies (e.g., 
Papovich et al. 2007; Rigby et al. 2008; Elbaz et al. 
2010; Nordon et al. 2010). 
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An additional uncertainty arises from the fact that a 
fraction of the hght emitted by these IR-bright galaxies 
may also be associated with embedded active galactic nu- 
clei (AGN). The separation of the AGN and star- forming 
component of a galaxy's IR output requires knowledge 
of its full SED, which is difficult to measure for a large 
number of high redshift systems. One technique to get 
at this separation of AGN and star formation activity is 
through the use of mid-infrared spectroscopy where one 
decomposes the spectra into aromatic feature and con- 
tinuum components (e.g., Sajina et al. 2007; Pope et al. 
2008; Murphy et al. 2009). For a heterogenous sample 
of IR-bright galaxies this technique has revealed that the 
contribution of AGN to the total IR luminosity output 
can span a large range, and is, on average, non- negligible 
being roughly ~30% (Murphy et al. 2009). 

In this paper we use deep 70 fim data from the Far- 
Infrared Deep Legacy (FIDEL; P.I.: M. Dickinson) sur- 
vey to provide improvcxl estimates for the; total IR lumi- 
nosities of 24 /im-selected sources in the Great Origins 
Observatories Deep Survey-North (GOODS-N; Dickin- 
son et al. 2003) field out to 2: » 2.8. Using these im- 
proved values, along with an empirical relation for the 
fractional AGN contribution to the total IR luminosi- 
ties of the 24 /im-selected sources, we show how the IR 
luminosity dcinsity evolves separately for AGN and star 
formation versus redshift. In doing so, we also look at 
how the contribution from populations of normal galax- 
ies {Lm < 10" Lq), LIRGs, and ULIRGs to the SFR 
density varies as a function of lookback time. 

The paper is organized as follows: In §2 we describe 
the sample properties and observations. IR luminosity 
estimates through SED fitting, along with a description 
of how we estimate contributions from AGN, are dis- 
cussed in §3. A look into the applicability of standard 
UV extinction estimates from the UV spectral slope for 
our sample of 24 /zm-selected sources is given in §4. A 
discussion of the results on the evolution of the IR lumi- 
nosity density versus redshift is given in §5, along with 
an explicit description of the uniqueness of this study. 
Finally, in §6, we summarize our conclusions. 

2. SAMPLE SELECTION AND 
MULTIWAVELENGTH PHOTOMETRY 

The GOODS-N field is located around the Hubble Deep 
Field-North at 12'^36'^55% -|-62°14'15" (J2000). In Fig- 
ure 1 we overlay source detections on the 70 ^m image 
of GOODS-N (Frayer et al. 2006) and provide a sum- 
mary of the source detections at 16, 24, and 70 /xm, as 
well as in the optical and hard-band (2.0 — 8.0 keV) X- 
rays in Table 1 (see Appendix A for details of existing 
radio data). The total number of sources detected in all 
3 infrared bands is 112. The multiwavelength data used 
here are the same that were used in the analysis of the 
mid-infrared spectroscopic sample presented in Murphy 
et al. (2009). 

2.1. Mid- and Far-infrared Spitzer Imaging 

Spitzer observations at 16, 24, and 70 /xm of the 
GOODS-N field were taken as part of various programs. 
Observations at 24 fim using the Multiband Imaging 
Photometer for Spitzer (MIPS; Rieke et al. 2004) were 
taken as part of the GOODS Legacy program (PI: Dick- 
inson) and reach an rms of ~5 fiJy (see Chary et al. 



TABLE 1 

Summary of Detections within ACS Coverage 



16 fj,m 24 /xm 70 /im Optical (zsso) Hard X-Ray 
Nags 1122 2196 133 1884 156 

Note. — Out of a total of 2664 24 ^um detected sources. 

2007). The catalog was created using IRAC prior po- 
sitions and sources having a 24 /xm flux density greater 
than 20 /iJy and a signal-to-noise ratio (S/N) greater 
than 3 were considered to be detections leading to a total 
number of 2664 detections; a 24 /im fiux density of 30 /iJy 
corresponds to a differential completeness of ~80% (Mag- 
nelli et. al. 2009). The 24 fim flux density among these 
sources spans a range between ~20— 4320 /iJy, more than 
a factor of ~200; the median flux density is ~70 /iJy. 
Calibration uncertainty at 24 fim is ^-^10%. 

Of the 2664 24 fim detections a total of 2196 sources 
are included within the optical {Hubble ACS; see §2.3) 
areal coverage of the GOODS-N field with 1111 spectro- 
scopic redshifts available (Cohen et al. 2000; Wirth et 
al. 2004; Barger, Cowie, & Wang 2008; D. Stern et 
al. 2010, in preparation). The differential spectroscopic 
completeness as a function of the 24 fim flux density is 
plotted in Figure 2, along with the corresponding 24 fim 
number counts. Among these sources having spectro- 
scopic redshifts, 1105 are at z < 3 having a median value 
of RsO.845 (see Table 2 for a breakdown using 6 redshift 
bins) . The 24 /im flux densities of these sources span a 
range from ~20 — 2440 fiJy (i.e., a factor of ~120), with 
a median value of ~80 fiJy 

Imaging at 16 fim was carried out using one of the 
two peak-up array cameras for the Infrared Spectrograph 
for Spitzer (IRS; Houck et al. 2004). Details of the 
imaging program and catalog can be found in Teplitz 
et al. (2005, 2010). The areal coverage of the 16 ^m 
map was slightly smaller than that of the GOODS-N 
24 fim map (see Figure 1). The rms noise of these data 
is ~6 fiJy and only sources having S/N ratios greater 
than 5 were considered to be firm detections leading to a 
total of 1297 sources detected. Of these, 1154 have 24 fim 
counterparts (1122 within the ACS coverage), while 775 
of these also had spectroscopic redshifts; 769 were at z < 
3. Among these 769 sources, their 16 fim flxix densities 
span a range from ~30 — 1420 /iJy (i.e., a factor of nearly 
~50) with a median value of 95 fiJy. The calibration 
uncertainly at 16 fim is ~10% and all non-detections 
were assigned the upper limit value of 30 // Jy. 

Deep MIPS 70 fim observations centered on, and ex- 
tending beyond, the GOODS-N fleld were performed as 
part of two separate imaging programs; GO-3325 (PI: 
Frayer et al. 2006) and the Far Infrared Deep Extra- 
galactic Legacy project (FIDEL; PI: Dickinson). The 
typical point-source noise of the 70 fim map is '-^^0.55 mJy 
(Frayer et al. 2006) which is roughly ^1.6 times larger 
than the confusion noise level of ctc = 0.35 ± 0.15 mJy 
(Frayer et al. 2009). Details of the data processing can 
be found in Frayer et al. (2006). The cataloged sources 
were considered detections if their 70 fim fiux density was 
> 2 mJy and they had a S/N > 6; these criteria corre- 
spond to a differential completeness of ~80% (Magnelli 
et. al. 2009) and a total of 167 sources are detected, 164 
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Fig. 1. — The location for each of the 2664 24 fim detected sources overlaid on the MIPS 70 /im map of the GOODS-N field. Corresponding 
detections at 16 ^im, 70 fira and 20 cm are shown; note that the areal coverage of the 16 /im map is smaller than the extent of the 24 /im 
observations. The blue outline shows the ft!l60 arcminute^ ACS areal coverage over which our analysis of the IR luminosity evolution 
versus redshift is conducted. 

were matched with an X-ray counterpart reported by 
Alexander et al. (2003). A total of 156 of these sources 
have firm detections in the hard band. 

The full band X-ray flux spans more than 3 or- 
ders of magnitude, ranging between 3 — 6950 x 
10~^^ erg cm~^ s~^, with a median flux of 87 x 
10~^^ erg cm~^ s~^. The hard-band X-ray detected 
sources span a smaller range in flux, spanning between 
9 — 4230 X 10"^'' erg cm~^ s^^, with a median flux of 
135 X 10^^'' erg cm^^ s^^. 

2.3. Optical and Near- Infrared Imaging 

The sal60 arcminute^ GOODS-N field has been im- 
aged extensively at optical wavelengths using the ACS 
camera on the Hubble Space Telescope in the following 
four filters: ^435, Veoe, ^7751 and zs5o (Giavalisco et 
al. 2004). In our analysis we use the catalogs for ver- 
sion 1 of the ACS data obtained from MAST. Additional 
ground-based imaging of GOODS-N in the J7-band was 



of which had 24 ^m counterparts (133 within the ACS 
coverage). Of the 164 sources matched with a 24 //m 
counterpart, 123 also had spectroscopic redshifts, all of 
which were at z < 3. The 70 /im flux densities range 
between ^2 — 19 mJy, spanning nearly an order of mag- 
nitude, and have a median value of 3.7 mJy. Sources 
not meeting the detection criteria were assigned the up- 
per limit value of 3 mJy. The calibration uncertainty at 
70 /im is ~20%. A summary of the total number source 
detections at 16, 24, and 70 /im within the ACS areal 
coverage can be found in Table 1. 

2.2. X-ray Imaging 

GOODS-N was observed using the Chandra X-ray 
Observatory for a 2 Ms exposure (Alexander et al. 
2003). The full band (0.5 - 8.0 keV) and hard band 
(2.0 — 8.0 keV) on-axis sensitivities are '^7.1 x 10^^^ and 
1.4xl0~^^erg cm~^ s~^, respectively. Of the 2196 24 fim 
detected sources included in the ACS areal coverage, 223 
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Fig. 2. — In the top panel the spectroscopic completeness for the 
24 /im sources included in the areal coverage of the ACS imaging 
is plotted against the 24 fim flux density. The completeness is 
differential, illustrating the fraction of sources with spectroscopic 
rcdshifts per each 24 fim flux density bin. The bottom panel shows 
total number counts of 24 /xm sources per flux density bin. 




Fig. 3. — A comparison between the photometrically recovered 
redshifts to the spectroscopic values. The distribution of relative 
errors dz = (^phot ~ ^spec) /(I + ^spec) has a median of Si —0.005 ± 
0.186. Outliers (i.e., those sources having \dz\ > 0.15) arc identified 
with open circles. Excluding these sources reduces the median and 
dispersion of the distribution of relative errors to fa 0.003 ±0.058. 



obtained using the prime-focus MOSAIC camera on the 
KPNO Mayall 4 m telescope. A subset of these data 
were described in Giavalisco et al. (2004) and Capak 
ct al. (2004) ; the present analysis uses a version of the 
[/-band data that is has approximately twice the expo- 
sure time. The ?7-band photometry was measured by 
detecting sources in the Subaru i?-band image of Ca- 
pak et al. (2004) and measuring {/-band fluxes through 
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total; 2065/2196 



ipec-z; nil 

phot-z: 865 
rac bump: 81 ^ 1 . 




Fig. 4. — Histogram of the rcdshift distribution for 24 ^m sources 
included in the ACS areal coverage. The maximum allowed pho- 
tometric redshift is set to 2phot = 3; such sources were not used 
in our analysis of the evolution of the infrared luminosity density 



matched apertures. Sources were considered detections 
in the UBViz bands if they had a S/N > 3, correspond- 
ing to a typical detection limit of 0.050, 0.035, 0.024, 
0.035, and 0.026 ^Jy, respectively. 

GOODS-N has also been observed at near-infrared 
(NIR) wavelengths with the Wide-Field Near Infrared 
Camera (WIRCAM) at the CFHT, including K,-band 
observing programs from Hawaiian and Canadian ob- 
serving programs, and J-band data from a Taiwanese 
program. WIRCAM has a field of view covering 21 f5 x 
21.'5, fully encompassing the GOODS-N ACS, IRAC and 
MIPS areas. Here we use reductions of the WIRCAM 
data from Lin et al. (2010, in preparation), consisting 
of 27.4 hours of integration time at J and 31.9 hours 
at Kg. The images have seeing with FWHM « 0'.'7. 
Sources were considered detections at J and Kg if they 
had a S/N > 3 and a flux density larger than 0.312 and 
0.456 /iJy, respectively. An independent analysis of the 
WIRCAM Kg-hanA data has recently been published by 
Wang et al. (2010). 

IRAC data at 3.6, 4.5, 5.8, and 8.0 /im were taken as 
part of the GOODS Spitzer Legacy program (M. Dick- 
inson et al. 2010, in preparation). The formal 3a limits 
for an isolated point source are 0.079, 0.137, 0.867, and 
0.951 /xJy for channels 1,2,3, and 4, respectively. How- 
ever, in practice, the signal to noise will depend on the 
degree of crowding with other nearby sources. 

2.4. Inclusion of Photometric Redshifts 

The 24 /im detections were cross-matched against the 
ACS BViz, CFHT JK^, and the 3.6 and 4.5 /im IRAC 
catalogs to provide photometric input for the Z-PEG 
code (Le Borgne & Rocca-Volmerange 2002) to com- 
pute photometric redshifts. Only sources having de- 
tections (i.e. S/N ratio > 3) in 4 or more of these 
wavebands are considered to have reliable photometric 
redshifts. A maximum redshift value of 3 was used 
while running Z-PEG. The distribution of relative errors 
dz = (2phot - 2spec)/(l + Zspec) among the 1050 sources 
having spectroscopic and photometric redshifts has a me- 
dian and dispersion of « —0.005 ±0.186 (Figure 3). Ap- 
plying an outlier criterion of \dz\ > 0.15, which excludes 
139 (ssl3%) of the sources, the median and dispersion of 
the relative errors are reduced to « 0.003 ± 0.058. 

Using all 1050 sources (i.e., including \dz\ > 0.15 out- 
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liers), we assign uncertainties to the photometric red- 
shifts by taking the standard deviation of the differ- 
ence between the spectroscopic and photometric redshifts 
((TAz) within each of the last 5 redshift bins listed in Ta- 
ble 2 resulting in uncertainties of 2; ± 0.33, 0.27, 0.28, 
0.46, and 0.46, respectively. These values, in addition to 
the imcertainties on the input flux densities in the SED 
fitting, were used in a Monte Carlo approach when es- 
timating uncertainties on IR luminosities derived with 
photometric redshifts (see §3.1.1). 

For an additional 81 sources not having enough broad- 
band detections to extract a reliable photometric redshift 
using Z-PEG; we use their peak in the IRAC channels to 
place them into a redshift bin based on the expected lo- 
cation of the 1.6 /xm stellar bump. Taking the updated 
Bruzual & Chariot (2003) stellar templates, which in- 
clude a revised prescription for stars on the asymptotic 
giant branch (AGB; see Bruzual 2007), extincted by 
an Ay = 1, it is found that the 1.6 /im stellar bump 
should shift into IRAC channels 2 and 3 for redshift bins 
of l.S < z < 1.9 and 1.9 < z < 2.8, respectively. For 
these sources, the center of the redshift bins are used 
when calculating IR luminosities (i.e., z « 1.6 and 2.35, 
respectively). Uncertainties in these redshifts were esti- 
mated by taking the standard deviation of the difference 
between the spectroscopic and IRAC-bascd photomet- 
ric redshifts where possible, resulting in uncertainties of 
.2 » 1.6 ± 0.47 and z w 2.35 ± 0.69. The choice of an 
Ay = 1 was for simplicity; we note that by increasing 
the assumed extinction to an Ay — 3, the center of the 
redshift bins decrease to z « 1.44 and 2.19, resulting in 
IR luminosities which are »i20% and 35% smaller, re- 
spectively, which is within our estimated errors. 

As with the IR luminosities derived using photometric 
redshifts above, uncertainties in the IRAC-based photo- 
metric redshifts were used to assign an additional uncer- 
tainty to the IR luminosity through a standard Monte 
Carlo technique. The redshift distribution for the 2065 
(out of 2196) 24 yum-selected sources included in the ACS 
areal coverage, for which we have a spectroscopic or pho- 
tometric redshift, is shown in Figure 4. The inclusion of 
these photometric redshifts results in a redshift complete- 
ness of «94%. We note here that photometric redshifts 
are only used for our analysis of the evolution of the IR 
luminosity density versus redshift. 

3. ESTIMATING STAR FORMATION AND AGN 
ACTIVITY 

In the following section we briefly describe our meth- 
ods for calculating infrared (IR; 8 — 1000 /im) luminosi- 
ties (Liji) for each 24 /im detected source. In addition, 
we describe a method to estimate the AGN contribu- 
tion to each IR luminosity. IR luminosities are used to 
estimate corresponding SFRs following the conversions 
given in Kennicutt (1998). The analysis presented here 
closely follows that of Murphy et al. (2009), where a 
more detailed description can be found; in Appendix A 
we present a similar comparison of IR, radio, and UV 
(1500 A) based SFRs, while in Appendix B we demon- 
strate that the bolomctric and AGN corrections given in 
Murphy et al. (2009) arc applicable to the present study 
of 24 /xm-selected galaxies. All calculations are made 
assuming a Hubble constant of 71 km s~^, and a stan- 
dard ACDM cosmology with Om = 0.27 and fl\ = 0.73. 



TABLE 2 
Sources per z Bin 



JV24 /im fim 



Spectroscopic Redshifts 



0.0 < z < 


0.4 


111 


30 


0.4 < z < 


0.7 


310 


43 


0.7 < z < 


1.0 


338 


28 


1.0 < z < 


1.3 


176 


12 


1.3 < z < 


1.9 


90 


6 


1.9 < z < 


2.8 


76 


3 


Spectroscopic and Photometric Redshifts 


0.0 < z < 


0.4 


123 


31 


0.4 < z < 


0.7 


333 


43 


0.7 < z < 


1.0 


430 


28 


1.0 < 2 < 


1.3 


333 


15 


1.3 < 2 < 


1.9 


349 


6 


1.9 < 2 < 


2.8 


362 


4 


Note. — 


Only 


sources 


within the ACS 



areal coverage axe considered. There axe 2 
sources with spectroscopic redshifts outside 
of the ACS areal coverage, as well as 4 ad- 
ditional sources with spectroscopic redshifts 
between 2.8 < z < 3.0. 

Throughout the paper average properties are measured 
by taking a median rather than a mean. 

3.1. Estimating IR luminosities of Galaxies 

In the following subsection we describe our methodol- 
ogy for deriving estimates of the IR luminosities among 
the sample. We also describe how we account for the 
fractional contribution of AGN to the IR luminosity of 
each source. 

3.1.1. Determining IR Luminosities from SED Fitting 

We fit the photometric data (i.e., the 16, 24, and 70 /zm 
flux densities where available) with the SED templates 
of Chary & Elbaz (2001) and then integrate between 
8- lOOO'/im. The SED libraries of Chary & Elbaz (2001) 
are used since they have recently been found to exhibit 
24/70 fxm flux density ratios which are better matched 
to the average observed for galaxies at 2; ~ 1 compared 
to the Dale & Hclou (2002) or Lagache, Dole, & Puget 
(2003) templates (Magnelli et. al. 2009). The best-fit 
SEDs are determined by a minimization procedure 
in which the amplitude (luminosity) and shape (tem- 
perature) of SED templates are allowed to scale freely. 
Fitting errors are estimated by a standard Monte Carlo 
approach using the photometric uncertainties of the in- 
put flux densities, as well as redshift uncertainties for 
sources having only photometric redshifts. In the cases 
where only upper limit values are available for the pho- 
tometric data, they are not incorporated into the formal 
minimization but are used to reject flts having asso- 
ciated flux densities greater than the upper limit. These 
IR luminosities are identified by LJr'^^''''^. 

For sources without far- infrared information (i.e., those 
sources not detected at 70 /im), which constitutes the 
majority of sources, ij^'^ ' values are estimated by av- 
eraging the IR luminosities from the best-fit Chary & 
Elbaz (2001) and Dale et al. (2001) SEDs. Since 
the Dale et al. (2001) SED templates are not nor- 
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Fig. 5. — In the left panel we plot the ratio of the 24 /im-derived IR luminosity to those estimated by SED template fitting with the 
16, 24, and 70 /im photometry for all 70 //m detections (open circles) versus redshift. Also included is the ratio of the 24 /im-derived to 
the best-fit IR luminosities for the mid-infrared spectroscopic sample of Murphy et al. (2009) (plus symbols). In the right panel we plot 
the same ratios against the 24 /Ltm-derived IR luminosity. ULIRGs (Lir > 10^^ Lq) are identified by asterisks. The vertical dotted line 
separates ULIRGs from non-ULIRGs and the dashed line is the ordinary least squares fit to the ULIRGs given by Equation 1. 



malized by luminosity, wc do so using a local trend be- 
tween IRAS colors and IR luminosity as the Dale et al. 
(2001) SEDs are described by a single parameter (i.e., 
/^(60 ^m)//y(100 A*m)) family. Due to the non-detection 
at wavelengths close to the peak of the far-infrared emis- 
sion, the IR luminosities of such sources may be consid- 
ered as upper limits. We note that for a heterogeneous 
sample (i.e., SMGs, AGNs, and optically faint sources) 
of 22 galaxies spanning the same redshift range being 
probed here, IR luminosities derived in this manner were 
found to be ^ 1.0 ±0.5 times those estimated from SED 
fitting with additional mid-infrared spectra and (for 18 of 
the 22 galaxies) 70 and/or 850 /im photometry (Murphy 
et al. 2009). 

We additionally calculate IR luminosities by fitting 
SEDs with the 24 iim photometry alone. This is also 
done by fitting both the Chary & Elbaz (2001) and Dale 
et al. (2001) SED template libraries independently, and 
averaging their integrated 8— 1000 /iin luminosities. The 
associated IR luminosity is designated as Lf^, and is, on 
average, within <5% of those derived from fitting either 
set of SED templates. The difference between the resul- 
tant Lj^ from the two sets of SED templates is charac- 
terized as a systematic error. 

3.1.2. Empirical Correction for 24 iim-derived IR 

Luminosities 

Murphy et al. (2009) have shown that IR luminosities 
derived by fitting local SEDs with 24 /xm photometry 
alone typically overestimate the true IR luminosity by 
an average factor of ^5 among sources having z > 1.4 
and Lj^ values > 3 x 10^^ Lq. We use that result to de- 
rive a empirical correction for Ljp^. In the left and right 
panels of Figure 5 we plot the ratio of 24- /im-derived to 
best-fit IR luminosities (ijp'^'''™) for all of the 70 /um 
detected sources having spectroscopic redshifts, versus 
redshift and 24 /um-derived IR luminosities, respectively. 
Also included in each panel is the ratio of 24-/xm-derived 
to best-fit IR luminosities for 18 galaxies presented by 
Murphy et al. (2009); the IR luminosities among these 
galaxies were fit using 70 and/or 850 /xm photometry 
and are considered fairly well determined. For the nine 
70 fim. detected sources which overlap between the IRS 



and larger GOODS-N samples, we use the IR luminosity 
reported by Murphy et al. (2009) because those esti- 
mates include the mid-infrarecl spectroscopic data in the 
fits. 

The trend reported by Murphy et al. (2009) ap- 
pears to hold. The 24 /xm-derived IR luminosities of 
galaxies having redshifts between 1.4 < z < 2.6 and 
■^iR ^ 3 X 10^^ Lq are overestimated by an average 
factor of ^4.6 ± 1.7. Galaxies at lower redshifts, hav- 
ing Lj]^ < 10^^ Lq, have 24-/im-derived IR luminosities 
which are consistent with our best-fit determinations us- 
ing additional longer wavelength data. These results are 
also consistent with recent Herschel findings confirming 
the reported overestimates in 24 /im-derived IR luminosi- 
ties by Murphy et al. (2009) (i.e., Elbaz et al. 2010; 
Nordon et al. 2010). 

While trends exist between the overestimates in the 
bolometric correction for both redshift and 24 /im- 
derived IR luminosity (see right panel of Figure 5), we 
believe the latter to be the dominant parameter as the 
trend between the ratio of 24 /zm-derived to best-fit IR 
luminosity versus 24 /im-derived IR luminosity appears 
tighter than versus redshift; the residual dispersion in 
the fit between the ratio of 24 /im-derived to best-fit IR 
luminosity and redshift (for z > 1.3) is ~40% larger than 
the fit versus luminosity, where if^ > 10^^ Lq . Further- 
more, individual lensed Lyman break galaxies (LBGs) at 
z ~ 2.5, having IR luminosities < 10^^ Lq, appear to 
show bolometric corrections that are similar to those of 
local galaxies (Siana et al. 2008; Wilson et al. 2008; 
Gonzalez et al. 2009). 

Consequently, we correct the 24 /im-derived IR lu- 
minosities using this information. For galaxies having 
< 10^^ Lq we do not apply a correction given that 
the median ratio of L^^ to the best-fit IR luminosity is 
1.1 ±0.27. For galaxies having Lf^ > 10^^ Lq, we apply 
an empirical correction by fitting these sources with an 
ordinary least square regression such that. 



STAR FORMATION AND ACCRETION HISTORIES AT HIGH REDSHIFT 



7 



TABLE 3 

Star Formation and Accretion Powered IR Luminosity Density Evolution 
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Note. — Vsfr(Mq yr-i Mpc^^) = 1.73 x lO"!" [pir(Lc5 Mpc-'-^) - pagnILq Mpc" 
may be considered as upper limits (see §3.1.3). *Perceiitagc of the SFR density. 
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Fig. 6. — The best-fit IR luminosities for each of the 1107 24 jum 
detections having spectroscopic redshifts of z < 3. Sources de- 
tected at 70 /im are identified by an open circle. 
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For the cases where these IR luminosity estimates exceed 
those derived using the additional 16 and 70 /xm data, we 
take the minimum value between these two values as the 
best estimate for the true IR luminosity of each source 
such that, 



LiR = min(L 



IR 



^IR 



(2) 



The fraction of cases for which Lj^'™'^'^ was much larger 

.1 r 16,24.70 /. t 24, corr / r 16.24,70 ^ n\ j r 

than Ljj^ (i.e., iyjp / J^m > 2) occurred tor 
~2% of the sample. Lastly, for the subsample of 22 
sources included in Murphy et al. (2009), we set iiR 
to their best-fit IR luminosities. 

The best-fit IR luminosities are plotted in Figure 6 as 
a function of rcdshift for all 1107 24 /im detected sources 
having spectroscopic redshifts of z < 3. The values range 
from 2.5 x 10^ — 1.1 x lO^'^ Lq, spanning more than 4 or- 
ders of magnitude, with a median of 7.2 x 10^" Lq. Lim- 
iting these to galaxies having firm 70 /ini detections, the 
median is increased by a factor of ~2.5 to 1.8 x 10^^ Lq 
and the range decreases to 1.0 x 10^ — 1.1 x 10^^ Lq. 
On the other hand, the range of the 24 /im derived 
luminosities, spans more than 5 orders of magnitude 
from 2.7 x 10* - 6.4 x 10^^ Lq, with a median value 
of 7.6 X 10^° Lq, roughly the same as the median of the 
best-fit IR luminosities. 



3.1.3. Constraints on AGN Fractions of IR Luminosities 

We make a statistical correction to the IR luminosities 
of each source by assigning a fractional contribution from 
AGN based on an empirical trend between AGN lumi- 
nosity and the difference between the 24 /zm-derived and 
true IR luminosities (Murphy et al. 2009) such that 



/ T AGN \ / r 24 \ 

log i-^—j = (0.73±0.08)log ( -^j+(2.5±1.0). (3) 



This relation relics on AGN luminosity estimates derived 
by first decomposing mid-infrared spectra into star for- 
mation (PAH template) and AGN (hot dust emission) 
components. Since hot dust associated with vigorous star 
formation may also contribute to the mid-infrared con- 
tinuum emission along with AGN, our estimates for the 
AGN luminosity could be considered upper limits. How- 
ever, this conclusion assumes that our choice of AGN 
template (i.e., Mrk 231) is appropriate for all sources, 
which may not be the case. Furthermore, it is cur- 
rently unclear as to whether the far-infrared luminosity 
of Mrk 231 is powered primarily by AGN (e.g., Armus 
et al. 2007; Veilleux et al. 2009; Fischer et al. 2010) or 
star formation (e.g., Downes & Solomon 1998; Farrah et 
al. 2003). If the latter is true, than that would again 
argue for the AGN IR luminosity estimates being upper 
limits. Correspondingly, the SFR density that we have 
derived would be increased by this amount. However, 
there is no reason to believe that these high-z sources do 
not contain AGN which contribute significantly to pow- 
ering the observed dust emission. Galaxies with SFR of 
~ 20 Mq yr~i occur in galaxies with stellar masses of 
~10i° Mq (Daddi ct al. 2007b) which should therefore 
harbor black holes of mass ^10^^* Mq if the black hole- 
stellar mass relation holds at z ^1 to 2 (McLure ct al. 
2006). Thus, our AGN estimates may be considered up- 
per limits, but in reality it is challenging to definitively 
quantify this without a precise knowledge of the AGN 
SED. 

By subtracting these values from our best estimates 
for the true IR luminosity, we can quantify the amount 
of IR luminosity arising from star formation alone such 
that Lf^f'^^ = Liu - Lf"^^. For the cases in which 
this empirical correction leads to AGN contributions that 
are larger than a source's estimated IR luminosity, we 
assume the_ source is completely powered by an AGN 



and set ifjp^ 
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3.2. The Evolution of the IR Luminosity Density 
versus Redshift 

Taking these improved estimates for the total IR lu- 
minosity of our 24 /im-selected sample of sources, we 
can investigate the evolution of the IR luminosity den- 
sity piK as a function of redshift. We do this in five dif- 
ferent redshift bins between 0.4 < z < 2.8 (i.e., the last 
5 redshift bins in Table 2). By using our IR luminos- 
ity estimates that have been corrected for the presence 
of extra infrared emission arising from embedded AGN, 
we convert the IR luminosity density into a SFR den- 
sity psFR to examine the evolution of the star formation 
history of the Universe over the last ~11.3 Gyr (see Fig- 
ure 7 and Table 3). In this analysis we only consider 
24 ^m sources included in the areal coverage of the ACS 
imagery since these were the sources targeted for spec- 
troscopic redshifts. The angular area of this region on 
the sky is approximately 160 arcminutes^. A detailed 
discussion of the evolution of the star formation and ac- 
cretion driven IR luminosity densities with redshift using 
our new results can be found in §5. 

3.2.1. Completeness Corrections for the Faint End of the 
Luminosity Function 

To correct for the additional luminosity arising from 
galaxies having 24 /im flux densities below the 24 /zm 
sensitivity limit (i.e., see Figure 6), we employ the lumi- 
nosity function derived from FIDEL data by Magnelli et. 
al. (2009, 2010). The luminosity function is parameter- 
ized by a broken power-law, whose evolution flattens once 
beyond z > 1.3, being consistent with a luminosity evo- 
lution proportional to {1 + z)^-^^^-'^ and (l + z)^'°^°'^ be- 
low and above z « 1.3, respectively. The corresponding 
density evolution for redshifts below and above z ~ 1.3 is 
proportional to {l+z)-°-^'^°-^ and (l-h^;)~^-^=^^-^, respec- 
tively. At high redshifts, these corrections will depend on 
the assumption that the infrared luminosity function has 
the same logarithmic slope that is measured locally (i.e., 
Sanders et al. 2003). 

By evolving the IR luminosity function appropriately 
for a given redshift, we integrate the luminosity contribu- 
tion for galaxies having IR luminosities between 10^ Lq 
and the IR luminosity corresponding to our 24 /im up- 
per limit of 30 /iJy at redshifts spaced by increments of 
Az = 0.05 within each of the 5 redshift bins between 
0.4 < z < 2.8. We then averaged the completeness cor- 
rections included in each of these redshift bins, weighted 
by the volume of each redshift increment. This value 
was then assigned as the effective completeness correc- 
tion for that redshift bin. The completeness corrections 
for these redshift bins are « 1.2, 1.8, 2.9. 3.8, and 
5.3 X 10® Lq Mpc~^, respectively. Similarly, we used 
this approach when accounting for corrections of galaxies 
in different luminosity classes (i.e., normal galaxies and 
LIRGs). The values of the completeness correction for 
LIRGs per each redshift bin are « 0.0, 0.0, 0.058, 0.70, 
and 2.6 x 10® Lq Mpc""^, respectively. The difference 
between the total completeness correction and these val- 
ues correspond to the completeness correction for normal 
galaxies alone. 

In Figure 8 we plot the space density evolution for 
LIRGs and ULIRGs. Also shown is the completeness 
corrected number density for LIRGs using the technique 



described above. These corrections are w 0.0, 0.0, 0.054, 
0.54, and 1.1 x 10~^ Mpc~^, respectively. 

3.2.2. Propagations of Errors for IR Luminosity Densities 

There are four primary sources of uncertainty in es- 
timating the IR luminosity densities: photometric red- 
shifts, photometric uncertainties, bolometric corrections, 
and completeness corrections. These are derived using a 
Monte Carlo technique from both photometry, and where 
applicable, the photometric redshift. The additional un- 
certainty of the scatter in our empirical correction for 
the overestimation of the 24 /im-derived IR luminosi- 
ties, which is fti27% for sources with Ljp < 10^^ Lq and 
w42% about the fitted regression line for sources with 
Lf^ > 10^^ Lq (see the right panel of Figure 5) is then 
added to the fitting errors in quadrature for each galaxy. 
Similarly, we include an uncertainty for the AGN sub- 
traction based on the residual dispersion (i.e., ^66%) in 
the empirical correlation (i.e., right panel of Figure 8 in 
Murphy ct al. 2009) used to derive Equation 3. This 
uncertainty is included in the quadrature sum above. 

The uncertainty in the completeness corrections for 
the faint end of the luminosity functions are based on 
the dispersion in the bolometric corrections (observed- 
frame Lir{/i^_L^(24 /im) ratios at each redshift bin be- 
tween 0.4 < z < 2.8 in Table 2) for the Chary & El- 
baz (2001) templates. Assuming effective dust temper- 
atures of Tdust = 36 ± 7 K, and a dust grain emissivity of 
/^dust = 1.6 (e.g., this temperature range has been esti- 
mated from average galaxy SEDs at these redshifts; Dye 
et al. 2007), we find that the corresponding bolometric 
correction uncertainties are 8, 16, 21, 32, and 62% for 
each redshift bin, respectively. These uncertainties are 
then added in quadrature to the uncertainties on the to- 
tal IR luminosity density of each redshift bin, which was 
based on the propagation of the uncertainties on the IR 
luminosities of each soiuce included in that redshift bin. 
The shaded regions plotted in Figure 7 illustrate the 1-a 
uncertainties. 

4. COMPARISON BETWEEN IR AND 
EXTINCTION-CORRECTED UV SFRs FOR 
24 /im-SELECTED SOURCES 

Having our improved estimates for IR luminosities, we 
now look to see how well local extinction corrections, 
based on the UV spectral slope, work for our sample of 
24 /zm-selected galaxies. A full treatment investigating 
local extinction correction among 24 /im-selected galax- 
ies in the BzK redshift range can be found in Appendix 
C. 

4.1. Typical Correction for UV Extinction 

UV-based SFRs (SFRuv) are calculated using the rest- 
frame 1500 A specific luminosities for each source using 
the U BViz photometry and the conversion given in Ken- 
nicutt (1998) as described in Murphy et al. (2009). We 
estimate the amount of extinction at 1500 A using the 
empirical relation between the UV spectral slope span- 
ning 1500 — 2600 A, defined as (3, to the total amount of 
extinction at 1600 A (Meurer et al. 1999), along with 
modeled extinction curves (i.e. Weingartner & Draine 
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Fig. 7. — The evolution of the comoving star formation and accretion-driven IR luminosity densities with redshift. In the top panel we 
plot pLi^ calculated using 24 based IR luminosities (asterisks) and our best-fit IR luminosities, fit using our other photometric data, 
after subtracting an estimate of the AGN contribution to the total IR luminosity (filled circles). An estimate of the AGN luminosity density 
is shown as crosses, which can be considered an upper limit. The total IR luminosity density (i.e., before subtracting out an estimate of 
the AGN contribution), corrected for sources below our sensitivity limits using the luminosity function work of Magnelli et. al. (2010), 
is given by the filled stars in each panel (see §3.2.1). For comparison with other work, we plot the IR luminosity density estimates of Le 
Floc'h et al. (2005, bold plus symbols) for our first 3 redshift bins; correcting our 24 /xm-derived IR luminosities with their IR luminosity 
functions (open circles), we find that our 24 /^m-derived IR luminosity density values agree, but are significantly above our best-fit SFR 
density determinations. In the bottom panel we show the total star formation driven IR luminosity density (i.e., after subtracting out an 
estimate of the AGN contribution), which have been corrected for sources below our sensitivity limits, as open stars. We also show the 
relative contributions to the star formation driven IR luminosity density (i.e., AGN corrected), calculated using all available photometric 
data, arising from low luminosity galaxies (diamonds), LIRGs (squares), and ULIRGs (triangles). The luminosity classi The 2 = data 
point for the total IR luminosity density, along with contributions from normal galaxies, LIRGs, and ULIRGs, was measured by integrating 
the z = luminosity function of (Magnelli et. al. 2009), while the 2 = AGN IR luminosity density was taken from Xu et al. (2001) and 
corrected for different cosmologies. The shaded regions include the la uncertainties. 
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Fig. 8. — The space density evolution of 24 ^m-selected LIRGs 
(squares), and ULIRGs (triangles). The contribution from LIRGs 
before applying a completeness correction (see §3.2.1) are given by 
the filled squares. Values at z = 0, along with the first two rcdshift 
bins for which we do not detect ULIRGs (i.e., z = 0.55 and 0.85) 
were taken from Magnelli et. al. (2009). 



OS 



10' r 

10^ r 
10' r 



0.66 



z < 2-6 



K 10" r 



10"' : 
10"^ : 

10"^ : 



P<®.-®..fi) ® 

(33 

If" 



O 70 detected 



0.0 



0.5 
/?(fi-I/) 



1.0 



Fig. 9. — We plot the ratio of IR derived SFRs estimated from 
SED fitting (plus UV derived SFRs uncorrected for extinction) to 
extinction corrected UV SFRs against the color excess E(B — V) 
inferred from fitting the rest-frame UV (1250 — 2600 A) slope us- 
ing UBViz optical imaging data; those galaxies which are 70 ^m 
detected are identified with circles. IR luminosities have not been 
corrected for AGN contributions, which makes little difference for 
the observed trend. 



2001; Draine 2003), such that, 



^1500 

mag 



4.65 -1- 2.09/3. 



(4) 



Since the redshifts of the sample galaxies are at z < 3, 
any correction to the UV slope due to absorption by the 
intergalactic medium is hkely neghgible. Sources having 
/3<— 3or/3>2 are excluded in this analysis as this 
is outside the range for which the UV spectra slope has 
been calibrated to extinction. These sources were gener- 
ally those having (low) redshifts for which /? was deter- 
mined using two data points that included the ground- 
based f/-band data, which is less sensitive and has much 
coarser resolution than the ACS imaging. The minimum 
redshift for which we have calculated a reliable extinc- 
tion correction UV luminosity is z w 0.66. Extinction 
corrected UV SFRs are denoted as SFR™^''. 

4.2. Applicability to 24 fJ-m Detected Sources 

In Figure 9 we plot the ratio of IR-I-UV to UV corrected 
SFRs versus the estimated E{B — V) color excesses (i.e.. 



E{B — V) = 0.123Ai5oo). For these comparisons with 
the UV SFRs, IR luminosities have not been corrected for 
the presence of AGN, however using the AGN subtracted 
values does not significantly affect the results. We note 
that the minimum redshift among these sources, where 
we could obtain a reliable estimate of the UV spectral 
slope, is z w 0.66. We have also limited this comparison 
to sources with z < 2.6. The UV corrected SFRs appear 
to be overestimated by more than a factor of >2, on 
average, among the 24 ^m-selected sources plotted, and 
we see a clear trend of decreasing ratio of IR+UV to 
UV corrected SFRs with increasing extinction suggesting 
that the UV slope is overestimating the extinction in a 
large number of sources. 

The fact that there is a trend at all suggests that the 
treatment for extinction correcting the UV flux densities 
is not adequate for our sample. To quantify this, we plot 
an estimate of ^1500 1 which is derived from the ratio of 
the total (IR + observed UV) to the observed UV SFR, 
against the UV spectral slope, /3 (see the left panel of 
Figure 10). A clear trend of increasing extinction with /3 
is found, for which an ordinary least squares fit results 
in the relation 



500 



mag 



2.75 ± 0.06 -f (0.75 ±0.04);3. (5) 



Our sample of 24 /im-selected galaxies typically fall well 
below the UV extinction correction of Meurer et al. 
(1999, i.e.. Equation 4), which is plotted as a dashed 
line. 

Taking these revised estimates for the UV extinction, 
we correct our UV luminosity and derive revised UV cor- 
rected SFRs (SFR^J'^-™"). In the right panel of Figure 
10 we plot the IR-I- (observed) UV SFRs against the re- 
vised, extinction corrected UV SFRs which exhibits a 
general one-to-one trend, albeit with a significant amount 
of (factor of ^3) scatter. Also shown in the right panel 
of Figure 10 is the result from fitting the IR-I- (observed) 
UV SFRs with the extinction corrected UV SFRs using 
the relation taken from Meurer et al. (1999). A clear dis- 
crepancy is found between these two fits suggesting that 
the extinction correction prescription given by Meurer et 
al. (1999) is not appropriate for our sample of 24 /xm- 
selected sources. Extinction corrected UV luminosities 
are a factor Ri2 larger using the local (starburst) relation 
of Meurer et al. (1999, i.e.. Equation 4) compared to our 
revised values using Equation 5. 

We believe this result may arise in part due to the 
age-extinction degeneracy among the lower redshift (i.e., 
0.66 ^ z < 1.4) objects in our sample. Attenuation rela- 
tions flatter than that of Meurer et al. (1999) have also 
been reported for samples of 'normal' star-forming galax- 
ies (e.g., Beh 2002; Kong et al. 2004; Buat et al. 2005), 
250 /im-selected galaxies at 2; < 1 (Buat et al. 2010), as 
well as UV-selected z ^ 1 galaxies when comparing to 
stacked X-ray SFRs (Laird et al. 2005). If the extinction 
is relatively lower in these objects, the UV slope will no 
longer be as sensitive to the amount of extinction due to 
the importance of variations arising from contributions 
by a galaxy's old stellar population. Another explanation 
for a flatter attenuation relation could be due to a steeper 
UV extinction curve among our sample of 24 /im-selected 
sources. While most normal z ~ 2 LBGs (i.e., UV- 
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Fk;. 10. In the left panel we plot our estimate of A1500, which is derived from the ratio of the total (IR + observed UV) to the observed 
UV SFR, against the measured UV spectral slope, ,3, for those galaxies shown in Figure 9. As in Figure 9, IR luminosities have not been 
corrected for AGN contributions. The over jjlottcd solid-line is the ordinary least squares fit to Aisoo versus /3 (Equation 5), while the 
dashed-linc indicates the UV extinction correction of Meurer et al. (1999, i.e., Equation 4). In the right panel we plot the total (IR + 
observed UV) SFR against our revised extinction corrected UV SFRs using the fit from the left panel. The dotted-linc is a one-to-one line. 
The dashed-linc indicates the fit to the IR -f- observed UV SFR versus the UV corrected SFR when using the UV spectral slope to correct 
for reddening given by Equation 4. This discrepancy clearly shows that the UV extinction correction using Meurer et al. (1999) is not 
appropriate for this sample of 24 /im-selccted sources. 



selected starbursts with 10^° L© < Lir < lO^^ Lq) fol- 
low the Meurer et al. (1999) relation, young (<100 Myr) 
LBGs tend to fall below the local starburst relation which 
has been attributed to a steeper, more SMC-likc, UV ex- 
tinction law in these sources (Reddy et al. 2006; Rcddy 
et al. 2010). This physical effect may also play a role in 
placing our sample of 24 ;um-selected galaxies below the 
Meurer et al. (1999) relation in the left panel of Figure 
10. 

5. DISCUSSION 

Following the results of Murphy ct al. (2009), we use 
deep imaging at 70 ^,m of GOODS-N, as part of the FI- 
DEL survey, to calculate the IR luminosities for 24 /xm 
detected sources. Using our improved bolometric cor- 
rections and estimates for the contribution of AGN (see 
§3.1), we characterize the evolution of the star formation 
and AGN powered IR luminosity density out to z ~ 2.8. 

5.1. Uniqueness of Sample and Analysis 

Here we highlight the major differences between our 
work for deriving the co-moving SFR and accretion en- 
ergetics over the past ^11 Gyr compared to others in the 
literature. First, we are working with a complete sample 
of 24 ;um-selected sources derived from the deepest 24 nm 
imaging carried out by Spitzer. The fact that these ob- 
servations were taken for the GOODS-N field provides 
a treasury of optical and NIR ancillary data for assign- 
ing redshifts (spectroscopic and photometric) to ?a94% or 
the sample. Additionally, deep 70 /im, and 850 /im data 
were available for a subset of the sample, allowing us to 
estimate the true IR luminosities of each source based on 
the deep 24 /im photometry alone (i.e., without having 
to rely on stacking analyses). Deep mid- infrared spec- 
troscopy has allowed us to empirically estimate the frac- 
tional AGN contribution to each IR luminosity, resulting 
in the first measure of obscured AGN energetics with red- 
shift, in addition to measuring contributions to the SFR 



density as a function of redshift for LIRGs and ULIRGs. 
Consequently, we have derived a new UV extinction cor- 
rection based on the measure of the UV spectral slope for 
24 jLtm-selected sources using these new, more accurate, 
estimates for IR luminosity. 

5.2. An Accounting of the Star Formation and 
Accretion Histories to z ~ 3 

In the top panel of Figure 7 we plot the IR luminosity 
density versus redshift for galaxies included in the last 5 
redshift bins of the bottom panel of Table 2 (i.e., 24 iim- 
selected galaxies covered in the ACS imaging). For the 
directly detected 24 fim sources, we plot independently 
the contribution from star formation (solid circles) and 
accretion processes (crosses) to our best-fit IR luminosi- 
ties. We also plot the IR luminosity density evolution 
one would predict based on tising the 24 /im-derivcd 
IR luminosities (asterisks). The corresponding IR lumi- 
nosity density estimates for our first three redshift bins 
taken from Le Floc'h et al. (2005) are also given (bold 
plus symbols); after applying a completeness correction 
based on their luminosity function work (open circles), 
we find that our 24 /im-dcrived IR luminosity density 
estimates are consistent. However, when the best-fit IR 
SFR is used (filled stars), we find that the Le Floc'h et al. 
(2005) results are significantly above our values. This is 
due to the more precise bolometric correction and AGN 
subtraction used in this paper. 

We estimate the total and star formation driven IR 
luminosity density evolution by including a complete- 
ness correction to our best-fit IR luminosities before 
(top panel; filled stars) and after (bottom panel; open 
stars) correcting for AGN (see §3.2.1); even with the 
completeness correction, the AGN-corrected values arc 
still well below what the 24 /im-derivcd IR luminosity 
density suggests. By using a more robust means of ex- 
cising out the AGN contribution to the total IR lumi- 
nosities, our value of the SFR density a,t z !v 2.35 is 
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PSFR « 0.21 ± 0.06 Mq yr~i Mpc"^. While this value 
is consistent with that reported by Reddy ct al. (2008) 
for a very similar redshift bin (i.e., 1.9 < z < 2.7), these 
authors did not account for a contribution from AGN. 
The fractional contribution of AGN to the total IR lumi- 
nosity density is found to decrease from <25% to <15% 
with increasing redshift between 0.55 ^ z < 2.35. On 
the other hand, the UV-corrected SFR density reported 
by Reddy et al. (2008), which assumes an average ex- 
tinction correction of a factor of 4.5, is « 12% larger than 
what is derived here. 

By averaging the SFR densities in our highest two red- 
shift bins, we obtain a z ^ 2 SFR density of psFR ^ 
0.17 ± 0.06 Mq yr^^ Mpc^'' which we can compare to 
other 2 ~ 2 values quoted in the literature. This value is 
slightly larger (i.e., «15%) than the z ^2 census value of 
0.15±0.03 Mq yr~^ Mpc^"^, computed from directly de- 
tected (optical, NIR, and submillimeter) galaxies (Reddy 
ct al. 2005), which can be considered a lower limit for 
the SFR density at this epoch. This value is also «50% 
larger than the z ^ 2 SFR density reported by Caputi 
et al. (2007), which is most likely due to their using a 
much flatter faint-end slope for their infrared luminos- 
ity function as pointed out by a number of authors (e.g., 
Reddy ct al. 2008; Rodighiero ct al. 2010; Magnelh et. 
ah 2010). 

In the bottom panel of Figure 7 we decompose the star 
formation driven IR luminosity densities into contribu- 
tions from normal galaxies (diamonds), LIRGs (squares), 
and ULIRGs (triangles). We note that sources arc first 
separated into IR luminosity classes (normal galaxies, 
LIRGs, and ULIRGs) based on their total IR luminosity 
(i.e., the IR luminosity before subtracting out an AGN 
contribution). By 2; f« 0.85 there is an equal contribu- 
tion to the SFR density from normal galaxies and LIRGs. 
The LIRGs continue to contribute equally with the pop- 
ulation of normal galaxies until z w 1.6, at which point 
the contribution of normal galaxies begins to decline out 
to our highest redshift bin at 2; « 2.35. We note that the 
completeness correction, which is unconstrained from the 
data presented here, is largest for normal galaxies in our 
highest redshift bin, reflected in the large uncertainties. 
The trend found for the space density evolution of LIRGs 
(Figure 8) is slightly different than their luminosity den- 
sity evolution. The space density of LIRGs appears to 
peak between 0.85 <z< 1.15 (« 1.6±0.1 x lO'^ Mpc'^) 
and remains flat until beyond z > 1.15 where it decreases 
to w 1.0 ± 0.4 X 10-3 Mpc-3 at z « 2.35. 

The ULIRGs, however, do not appear to account for 
more than ~20% of the SFR density until redshifts of 
z >1.6, at which point they become roughly comparable 
with the contribution from LIRGs once beyond z > 2, 
contributing at the «42% level by z « 2.35 (sec Table 
3). Our result is larger than the 30% contribution of 
ULIRGs to the total SFR density reported by Reddy et 
al. (2008) in this redshift range (i.e., 1.9 < z < 2.8). 
On the other hand, the total contribution from LIRGs 
(37%) and ULIRGs (33%) at z ~ 2 is «70%, signifl- 
cantly smaller than the w93% contribution from LIRGs 
(48%) and ULIRGs (45%) reported by Rodighiero ct al. 
(2010). Similarly, our result is in contrast with claims 
that ULIRGs alone dominate the SFR density at z ~ 2 
(e.g., Perez-Gonzalez et al. 2005). We attribute these dif- 



ferences between previous studies to our improved IR lu- 
minosity estimates and accounting for AGN, which pref- 
erentially affects the most luminous 24 /xm sources. In 
the case of Rodighiero et al. (2010), these authors re- 
moved type-1 quasars in their SFR density calculations, 
assuming the luminosities of such sources are entirely ac- 
cretion driven, but did not include a contribution from 
the bulk of the population as done here: they note that 
the removal of these sources did not have a large im- 
pact on their reported SFR densities, and is thus not the 
reason for differences between our results. 

Similar to Magnelli et. al. (2010), we find that the 
IR luminosity density of LIRGs remains roughly flat be- 
tween 1.15 < z < 2.35, however we also find the con- 
tribution from ULIRGs, and the total SFR density, may 
increase with redshift over this range. These authors also 
find that ULIRGs do not dominate the SFR density at 
z ~ 2, albeit they report an even smaller contribution 
from ULIRGs (i.e, 17%) than the 33%. reported here. 
While ULIRGs appear to contribute significantly to the 
SFR density at z w 2.35, we note that the contribution 
made by ULIRGs to the total SFR density of the Uni- 
verse integrated between < z < 2.35 is small, being 
«20%, where as LIRGs and normal galaxies contribute 
equally at the Ri40%i level. Furthermore, as shown in Fig- 
ure 8, the space density of ULIRGs is significantly less 
(i.e., a factor of ps3) than LIRGs at z « 2.35. 

While the contribution of ULIRGs appears to be large 
in our final redshift bin, it is not clear that this trend 
should continue to higher redshifts. Recent work by 
Magdis et al. (2010) report on the IR luminosity 
density of UV-selected ULIRGs at z ~ 3, being ~ 
7.5ttl ^ ^0 Mpc"^. This number is significantly 
smaller than (i.e., «15% of) our inferred SFR density 
for z 2.35 ULIRGs, and provides a lower limit on the 
fractional contribution of ULIRGs to the SFR density at 
z ~ 3 since the UV selection is likely to miss more ob- 
scured, UV-faint ULIRGs. We also note that the space 
density of the 24 /im-selected ULIRGs at z « 2.35 (i.e., 
w 3.0 ±0.6 X 10"'' Mpc-3) jg a fj^pi-oj. of ^20 times larger 
than the space density of z ^ 3 UV-selected ULIRGs 
reported by Magdis ct al. (2010). 

It is possible that our final redshift bin may be over 
predicting the SFR density due to an vmder estimation 
of AGN at these epochs since we are assigning the AGN 
contribution based on a relation which was derived for a 
heterogenous sample of only Ri20 galaxies (Murphy et al. 
2009). However, even by assuming that the IR luminosi- 
ties estimated for all hard X-ray detected sources, along 
with those sources having an IRAC SED described by a 
power law, arc completely AGN powered, this does not 
flatten the SFR density evolution nor bring the ULIRG 
contribution below that of the normal galaxies in our 
highest redshift bin; the SFR density for the flnal redshift 
bin decreases to psFR ~ 0.19 Mq yv~^ Mpc-^ while the 
fractional contribution by ULIRGs decreases to ~39%. 

While our results suggest an increase to the IR lumi- 
nosity density out to z«2.35 and that ULIRGs may con- 
tribute significantly to the SFR density beyond z > 2, 
we note that the uncertainties in our final redshift bin 
are quite large. Deep surveys with Herschel, such as the 
GOODS- Herschel open time key project (PI: D. Elbaz) 
will likely settle this ambiguity by properly measuring 
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the peak of the FIR SEDs for LIRGs and ULIRGs out 
to redshifts of 2; ~ 2 and 2 ^ 4, respectively. 

6. CONCLUSIONS 

In the present study wc have built on the findings of 
Murphy ct al. (2009) to determine how the star forma- 
tion and aeeretion driven IR luminosity densities have 
evolved over the; last «11.3 Gyr. This has been made 
possible through improved estimates of IR luminosities, 
achievable by further constraining fits to local SED tem- 
plates using deep MIPS imaging at 70 /xm taken as part 
of FIDEL. Our conclusions can be summarized as follows: 

1. IR (8—1000 iim.) luminosities derived by SED tem- 
plate fitting using observed 24 /um flux densities 
alone overestimate the IR luminosity by a factor 
of ^4 when their 24 /im-dcrived IR luminosity is 
> 3 X 10^^ Lq. This discrepancy appears to be the 
result of high luminosity sources at z ^ having 
far- to mid-infrared ratios, as well as aromatic fea- 
ture equivalent widths, which are more typical of 
lower luminosity galaxies in the local Universe. 

2. After accounting for star formation and AGN con- 
tributions to the total IR luminosities, we find that 
AGN and star formation activity appear to roughly 
track one another with AGN typically accounting 
for <18% of the total IR luminosity density inte- 
grated between < z < 2.35. The AGN fraction 
slightly decreases from <25% to <15% with in- 
creasing redshift between 0.55 ^ z < 2.35. Our ob- 
servations also hint that the star formation driven 
IR luminosity density may increase with redshift 
between 1.15 ^ z < 2.35, however, within uncer- 
tainties, our results are consistent with a flat evo- 
lution over this redshift range. 

3. The SFR density (i.e., IR luminosity density cor- 
rected for AGN contamination) is dominated by 



normal galaxies and LIRGs at comparable levels 
(i.e., each at « 40 - 50%) between 0.85 <z< 1.6. 
LIRGs continue to contribute at a similar level out 
to z Ki 2.35 which is in contrast with the ULIRGs, 
which transition to becoming a signiflcant contrib- 
utor to the SFR density (i.e., comparable with 
LIRGs) only once beyond z >2. 

4. Among our sample of 24 /xm-selected sources, we 
find that local prescriptions used to estimate UV 
extinction corrections based on the UV spectral 
slope typically overestimate the true extinction by 
a factor of >2. Accordingly, using our improved 
estimates for IR luminosity, we have derived a new 
UV extinction correction based on the UV spectral 
slope which may be more appropriate for 24 /^m- 
selected sources. 
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APPENDIX 

A. IR LUMINOSITY AND SFR COMPARISONS 

Miirphy et al. (2009) compared IR luminosities derived from 24 jim data alone to those derived when additional 
measurements are available, such as 16 and 70 /im photometry, submillimeter data, and mid-infrared (IRS) spec- 
troscopy. Here, we extend this to a much larger sample of GOODS-N sources out to z < 3. The redshift bins used, 
along with the number of sources per bin, per selection criteria, are given in Table 2. Only sources having spectroscopic 
redshifts are considered. Through this comparison, discrepancies between 24 /um-derived and best-fit IR luminosities 
are investigated. Although the 70/im data are less sensitive than the 24 /xm data, they help identify biases in the 
bolometric corrections at the bright end of the luminosity function. 

A.l. Trends with Redshift 

In Table 4 we give the median ratio of the 24 /im-derivcd, 24 /im-corrected, and radio-derived IR luminosities 
(see §A.3.1) to the best-fit IR luminosities per redshift bin along with measured dispersions. The 24 /im-derived IR 
luminosities appear to do well by matching our best-fit IR luminosities up to redshifts of z < 1.3, at which point the 
24/im-derived IR luminosities appear to overestimate the true IR luminosity by > 30%. More interestingly, we find 
that the dispersion in the ratios jumps significantly once past a redshift of 2; > 1.3, increasing by a factor of ~2 and 3 
in the last two redshift bins, respectively. 

If we only consider those sources detected at 70 (i.e., the brightest sources) we flnd that the discrepancy between 
the 24 /im-derived and best-fit IR luminosities is even larger, along with the measured dispersions, once beyond z > 1.3. 
The 24 /im-derived IR luminosity is a factor of ^ 2 and 6 times larger than the best-fit IR luminosities in the two 
highest redshift bins. These discrepancies are similar to what was reported by Murphy et al. (2009), whose sample 
was selected to have /i/(24 /xm) > 200 /xJy and mid-infrared spectroscopy. This result is also consistent with that of 
(Papovich et al. 2007) whom found 24 fim derived IR luminosities to be a factor of 2 to 10 times larger than those 



14 



MURPHY ET AL. 



TABLE 4 

BoLOMBTRic Correction Adjustments per z Bin 
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Note. — Only sources having spectroscopic redshifts considered. 



TABLE 5 

Bolometric Correction Adjustments per 24 //m-Derived IR Luminosity Bin 
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Note. — Only sources having spectroscopic redshifts considered. 

derived from stacking 70 and 160 fiia data for a sample of 24 /im bright (i.e., /i/(24 /im) > 250 /iJy) galaxies lying 
in a redsliift range between 1.5 < z < 2.5. Consequently, it appears that luminosity, and not redshift, may be the 
more important parameter associated -with the improper SED fitting -when using 24 /xm photometry alone. This result 
agrees with what was shown in §3.1.2 where the residual dispersion between the ratio of the 24 /^m-derived and best-fit 
IR luminosities was smaller when fitting versus luminosity as opposed to redshift. 

In comparing the radio-derived to best-fit IR luminosities among the 20 cm detected sources, we find that the radio 
(via the FIR-radio correlation) does not yield reliable IR luminosity determinations as they are >2 times larger, on 
average, in each redshift bin. This result is likely due to the difference in the depths of the radio and Spitzer surveys 
with the radio survey being much shallower, and thus less sensitive to galaxies of the same luminosity compared with 
the 24 fxm data (see §A.3.2). 



A. 2. Trends with 24 nm-Derived IR Luminosities 

In Table 5 the same IR luminosity ratios are given as in Table 4, but per 24 /xm-derived IR luminosity bin instead of 
redshift. For values of L^^ > lO^^.S the 24 /um-derived luminosity is a factor of ~4 times larger than the best-fit 

values, on average. Not too surprisingly, when only considering those sources detected at 70 /im, we find nearly the 
same trend. A comparison between the radio-derived and best-fit IR luminosities among the 20 cm detected sources 
shows that for all luminosity bins where Ljp > 10^° Lq, the radio-derived luminosities are >2 times larger than the 
best-fit estimates. Interestingly, we find no trend with increasing luminosity bin. 

A. 3. Comparison with Radio-Derived SFRs: The FIR-Radio Correlation versus z 

The GOODS-N field has been imaged at 1.4 GHz using the Very Large Array for a total of 165 hours in all 4 
configurations (Morrison et al. 2010). Since the VLA is centrally condensed, the observing time per configuration was 
scaled as follows, A-array (1), B-array (1/4), C-array (1/16), and D-array (1/64) (Owen & Morrison 2008). Using 
such an integration scaling provides the best sensitivity for extended sources. The final radio mosaic has a local RMS 
of ^3.9 /iJy near the phase center covering the GOODS-N ACS area. 

Taking the published GOODS-N 1.4 GHz catalog of Morrison et al. (2010), where 1230 discrete sources have been 
detected above a 5(t threshold, a total of 342 of lie within the ACS+24 /xm coverage. Of these sources, 286 could be 
matched to 24 /xm detections. Among these, 194 had spectroscopic redshifts, and 192 were at z < 3. The 20 cm flux 
densities among these 192 sources span a factor of ~33, ranging from ~21.3 — 704 /iJy. The median 20 cm flux density 
is ~46 /uJy. Sources not matched with a radio counterpart were assigned the upper limit value of 19.5 /xJy. 

A.3.1. Radio- Derived SFRs 

A nearly ubiquitous correlation is known to exist between the far-infrared (FIR; 42.5 — 122.5 /um) dust emission and 
predominantly non-thermal (e.g., 1.4 GHz) radio continuum emission arising from star- forming galaxies (de Jong et 
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TABLE 6 
IR-Radio Correlation per z Bin 
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Note. — Only sources having spectroscopic redshifts considered. The logarithmic 
IR/radio ratios and associated dispersions are given in units of dex. 



al. 1985; Helou et al. 1985). The most prominent feature of this correlation is that, at least in the local Universe, it 
extends for galaxies spanning nearly 5 orders of magnitude in luminosity while exhibiting a scatter which is less than 
a factor of ^2 (e.g., Yun, Roddy, & Condon 2001). Massive star formation provides the common link relating these 
two emission processes both globally and on <kpc scales within galaxy disks (e.g., Murphy et al. 2006; Hughes et 
al. 2006; Murphy et al. 2008). Consequently, optically thin radio continuum emission is often considered to be a 
very good SFR diagnostic. While the FIR-radio correlation has largely been established in the local Universe, there 
are hints that it may remain constant out to high redshifts (e.g., Garrett 2002; Appleton, et al. 2004; Frayer et al. 
2006; Kovacs et al. 2006; Beelen et al. 2006; Vlahakis et al. 2007; Murphy et al. 2009; Murphy 2009; Sargent et 
al. 2010a,b; Bourne et al. 2010). 

We parameterize the IR-radio correlation following a similar quantitative treatment of the FIR-radio correlation (i.e. 
Helou et al. 1985), except that we use the total IR (8 — 1000 jim) luminosity, rather than the FIR fraction such that, 

( 3.75x101^.(20 cm) )- (^1) 
Rest-frame 1.4 GHz radio luminosities are calculated such that, 

i.(20 cm) = A'kDIS^{2Q cm)(l + z)"''^ , (A2) 

which includes a bandwidth compression term of (1 + z)~^ and a K-correction of (1 + z)°' to rest frame 1.4 GHz. This 

assumes a synchrotron power law of the form Su oc v~'^ with a spectral index a for which a value of ^0.8 is assumed 
(Condon 1992). The median (Jir value reported for a sample of 164 galaxies without signs of AGN activity is 2.64 
with a dispersion of 0.26 dex (Bell 2003). 

By rewriting Equation Al, we can use the radio measurements to independently obtain estimates of the IR luminosi- 
ties of each source for comparison with our template fitting. For q = 2.64 the radio-based IR luminosities are defined 
as 

Lf^ = 1.64 X 10^^L.(20 cm) (A3) 

We assign a factor of ~2 uncertainty to these luminosities since this is the intrinsic scatter among (e.g., Yun, Reddy, 
& Condon 2001; Bell 2003) and within (e.g.. Murphy et al. 2006, 2008) star-forming systems in the local Universe. 

A. 3. 2. Lack of Evolution with Redshift out to z <^ 3 

In Table 6 we give the median IR/radio ratio for each subset of sources per redshift bin given in the top of Table 
2 along with the associated dispersion. By simply taking all 1107 24 /xm detections with spectroscopic redshifts less 
than 3 we find median IR/radio ratios which are typically more than a factor of 3 lower than the locally measured 
value with a dispersion > 2 for each redshift bin. This is likely driven by the fact that only 193 (i.e., ~17%) of these 
sources are radio detected while the remainder have IR/radio ratios calculated using the radio upper limits and are 
intrinsically much fainter. 

Focussing on the subset of sources for which we have firm detections at 70 /xm, and are most confident in the 

determination of their estimated IR luminosities, we find that their IR/radio ratios also show no clear signs of evolution 
with redshift while displaying a dispersion that is a factor of <2 for each redshift bin. In looking at the lowest redshift 
bin for these sources, we find the median IR/radio ratio and dispersion to be 2.62 and 0.28 dex, nearly identical 
to the average value of 2.64 dex reported by Bell (2003). At the highest redshift bin, the median qiR value also 
matches that for what measured in the local Universe. While the dispersion in qi^ among all 70 /Ltm detected sources 
is found to be as small as what is measured in the local Universe, being ~0.24 dex, we do note that the median value 
{qm, = 2.50 dex) is a factor of ~1.37 (0.14 dex) lower than the value reported by Bell (2003). Since this departure is 
much smaller than the la dispersion in the locally measured IR-radio correlation, and smaller in magnitude compared 
to the errors associated with the SED template fitting, speculating on possible physical mechanisms driving this offset 
is not warranted. 

Instead, taking the IR/radio ratios for only sources having 20 cm detections, we find that the median ratios are a 
factor of ~ 1.7 — 3 times lower than what is measured in the local Universe. The dispersion in for these sources 
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Fig. 11. — The ratio of 8 to 24 /im flux densities as a function of redshift for all 24 fim sources detected at 8 ^m, and for which a redshift 
(spectroscopic or photometric) could be derived. Bold filled circles indicate the 22 galaxies included in the IRS study of Murphy et al. 
(2009); of these 22 galaxies, those which are SMGs and detected in the hard-band (2.0 — 8.0 keV) X-rays are identified by a star and cross, 
respectively. 



is also generally >2 in each bin, and increasing with redshift. As with the radio derived SFRs for all 20 cm sources, 
wc attribute this systematic deviation in the IR/radio ratios to arise from the fact that the 20 cm detections are 
the brightest objects at each redshift for which the radio emission may have a considerable contribution from an 
AGN leading to overestimates in SFRs. Recent work by Sargent et al. (2010a) has shown that this apparent negative 
evolution in the FIR-radio correlation arises from comparing IR and radio selected samples having dramatically different 
sensitivity limits, with IR observations typically being much deeper than data from current radio surveys. Taking this 
into account, Sargent et al. (2010a) have used a survival analysis to demonstrate that the FIR-radio correlation out 
to similar redshift is consistent with the canonical value. Thus, we do not see any clear signatures of evolution in the 
FIR-radio correlation with redshift among the 70 /um detected sources in GOODS-N, nor any significant changes in 
the scatter per redshift bin, although the number of sources having spectroscopic redshifts for 2: > 1.3 is highly limited. 
This is not surprising since any expected evolution with redshift due to increased energy losses to cosmic-ray electrons 
from inverse Compton scattering off of the CMB, resulting in depressed synchrotron emission from galaxies, should 
not become important until beyond z > 5 (e.g., Murphy 2009; Lacki & Thompson 2010) 

B. APPLICABILITY OF IRS SAMPLE RESULTS TO THE FULL 24 ^m SAMPLE 

We compare the mid-infrarecl photometric properties of the 22 galaxies studied by Murphy et al. (2009) to the 
entire 24 /im-selected sample. This is done to demonstrate that applying the bolometric and AGN corrections derived 
in Murphy et al. (2009) to the present sample of 24 /im sources is justified. In Figure 11 we plot the observed ratio 
of 8 to 24 /zm flux densities versus redshift for each 24 ^m source detected at 8 /xm, and for which a spectroscopic 
of photometric redshift could be measured (i.e., 2034 of 2196 24 fim detected sources in the ACS coverage). The 22 
galaxies from Murphy et al. (2009) are shown using bold filled circles. Sources which are SMGs, and presumably star 
formation dominated, are identified by stars while sources which are hard X-ray detected and likely harbor an AGN, 
are identified by crosses. We find that the sample of Murphy et al. (2009) spans nearly all of the 8/24 /xm phase 
space covered by the full 24 fxm selected sample independent of the fact that the IRS sample was flux limited down to 
200 iiJy. 

A common diagnostic tool for identifying AGN is the IRAC color-color plot (e.g.. Lacy et al . 2004; Stern et al. 
2005). While commonly used as a means to identify AGN, it is worth pointing out that this diagnostic does have 
its shortcomingings. Contamination of the AGN "wedge" will occur by both low- and high-z star-forming systems, 
which have similar IRAC colors (e.g., Barmby et al. 2006; Donley et al. 2007; Cardamone et al. 2008). In Figure 12 
we create such a plot for all IRAC detected 24 /im sources (i.e.. 2396 of 2664 24 fim detected sources) to see if the 
IRS sample only occupies a single region of this phase space. We again identify the sample of Murphy et al. (2009) 
using bold filled circles, while SMGs and hard-band (2.0 — 8.0 keV) detected sources are plotted using a star and 
cross, respectively. It is clearly shown that the IRS sample is well distributed in IRAC color-color space, touching 
locations occupied by AGN and star-formation dominated systems. Consequently, it appears that the IRS sample 
is representative of the entire 24 /xm population, suggesting that applying the improved bolometric corrections and 
estimates for the fractional IR luminosity contributions from AGN given in Murphy et al. (2009) is warranted. 

C. SFR COMPARISONS AMONG BzK SOURCES: OCCURRENCES OF IR EXCESS SOURCES 

Daddi et al. (2007a) have identified discrepancies between the UV, radio, and, IR derived SFR estimates in z ~ 2 
galaxies selected using the BzK selection technique. Following the criterion of Daddi et al. (2007a, b), we calculate 
whether galaxies have a mid-infrared excess by measuring the ratio of IR-|- (observed) UV SFRs to that of the extinction 
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Fig. 12. — An IRAC color-color plot for all IRAC detected 24 /^m sources usually used to identify galaxies hosting AGN (e.g., Stern et 
al. 2005). Bold filled circles indicate the 22 galaxies included in the IRS study of Murphy et al. (2009); of these 22 galaxies, those which 
are SMGs and detected in the hard-band (2.0 — 8.0 keV) X-rays are identified by a star and cross, respectively. The distribution of the 
IRS sample fills in most of the phase spax;e covered in the IRAC color-color plot, indicating that the sample appears quite representative 
of galaxies spanning a range of types. 

corrected UV SFR. Sources having 

/ SFRiR + SFRuv\ ^ n c //-ii \ 

[ SFR-.. j>0.5 (CI) 

arc considered to be "mid-infrared excess" sources. We note that while the numerator contains a term for the UV 
emission which may escape the galaxy before being absorbed and re-radiated in the infrared by dust, this term is often 
negligible compared to the IR-based SFR term. 

Using the 24 /zm-derived, radio-derived, and best-fit estimates for IR luminosity to derive SFRs, we plot the ratio of 
IR-I-UV to UV corrected SFRs for sources having spectroscopic redshifts in the range between 1.4 < z < 2.6 detected 
at 24 /im for which we were also able to estimate a rest-frame 1500 A flux density. The IR based SFRs in the top panel 
of Figure 11 are calculated using the 24 /im-dorived IR luminosities, the middle panel uses our best- fit IR luminosities, 
and the bottom panel uses IR-based SFRs estimated using the radio continuum imaging and the FIR-radio correlation. 

Excluding the hard band (2.0 — 8.0 keV) X-ray detected sources, for which AGN are likely present, we find that the 
number of sources lying above the mid- infrared excess criterion decreases dramatically (i.e., from 7 to 4, a factor of ^-^2) 
by using the best-fit IR luminosities rather than those derived using the 24 /im photometry alone. This is consistent 
with the findings of Murphy et al. (2009) who reported that Ri50% of their "mid-infrared excess" sources could be 
accounted for by better constrained bolomctric corrections rather than the subtraction of emission from obscured AGN 
which was estimated using deep mid-infrared spectroscopy. Instead, looking at the radio derived IR luminosities for 
those sources detected at 20 cm, we find that the number of mid-infrared excess sources is actually larger than when 
using the 24 fim derived IR luminosities, being 9 and 6, respectively. 

C.l. Corrections for AGN 

While the improved bolometric correction was able to account for f«50% mid-infrared excess sources (Figure 13), we 

find that by subtracting our estimates for the AGN contribution to the IR luminosities, all galaxies remain mid-infrared 
excess sources (top panel of Figure 14). The AGN contribution is indicated by the distance between the horizontal 
line and solid points. This finding is consistent with Murphy ct al. (2009) who reported that the existence of the 
mid-infrared excess sources was dominated by an overestimate in the bolometric corrections rather than embedded 
AGN. This is contrary to the suggestion inferred by Daddi et al. (2007a), whose conclusion was based on the presence 
of a hard X-ray source through X-ray stacking. Accordingly, our result suggests that the sky and space densities of 
Compton-thick AGN reported by Daddi et al. (2007b) are high by a factor of ~2, and are likely more consistent with 
values of ~1600 deg^^ and ^1.3 x 10^** Mpc^'', respectively. 

In the bottom panel of Figure 14 we plot the ratio of AGN corrected IR (plus the observed UV) SFRs to extinction 
corrected UV SFRs and find the same occurrence of mid-infrared excess sources (4 sources, 3 of which are 70/im 
detected), again suggesting that AGN are not causing the discrepancy between the IR and UV SFRs; among all the 
galaxies plotted, the application of Equation 3 suggests that AGN only contribute <25% of their total IR luminosity, 
on average. To explain the galaxies which persist as mid-infrared excess sources. Murphy et al. (2009) suggested 
that it may have to due with improper extinction corrections relying on the slope of the UV continuum such that the 
extinctions for these sources are underestimated, thereby yielding underestimates for the UV-derived SFRs. In this 
scenario the extinction within a galaxy is so high that the ISM becomes optically thick, causing the relation between 
extinction and the slope of UV continuum to fail. 
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Flo. 13. Wc plot the ratio of IR- based SFR estimates (plus the observed UV derived SFR) to extinction corrected UV SFRs as a 
function of redsiiift between 1.4 < z < 2.6 for sources which wc were able to measure rest-frame 1500 A flux densities. In the top panel 
the IR based SFR. was calculated using IR luminosity estimates by SED fitting the 24 /xm flux densities alone while the middle panel uses 
the best-fit IR luminosities to calculate the IR-based SFRs. In the bottom panel we use the 20 cm radio continuum flux densities and the 
FIR-radio correlation to derive the SFRs (SFR^^^:). The SFRs are considered to be upper limits for those galaxies not detected at 20 cm 
(downward arrows). In each panel the dashed line indicates where the logarithm of this ratios is equal to 0.5 dex; galaxies having values 
higher than this are considered to be IR-excess sources (Daddi et al. 2007a). Galaxies detected at 70 //m axe identified with open circles 
while those that are detected in hard band (2.0 — 8.0 keV) X-rays are indicated by a cross. 



C.2. UV Extinction Corrections 

To compare how well simple color relations work to estimate extinction versus actually fitting the UV slopes of 
each objects, we plot the ratio of IR-f-UV to UV corrected SFRs in the in Figure 15 against E{B — V) color excesses 
calculated using the color relation given by Daddi et al. (2004, E{B — V)jjo4 = 0.25(B — ^-|-0.1)^s; left panel) and from 
a proper fit to the rest-frame UV continuum (right panel). The UV corrected SFRs used in the ratios were calculated 
using the respective methods to compute the extinction corrections. Figure 16 also compares the color excess E{B — V) 
derived by the two methods. The direct fitting method tends to yield larger reddening values, particularly for objects 
with redder UV spectral slopes. The median E{B — V) changes from ~0.25 to 0.29 mag between the color-fitting BzK 
method of Daddi et al. (2004) and the direct fits to the UV slope for individual galaxies as defined here. The median 
UV extinction correspondingly changes from ~2.0 to 2.4 mag, corresponding to an increase in the median UV-derived 
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Fig. 14. — The same as the top two panels in Figure 11 except that wc correct the 24 /im and best-fit IR luminosities for the contribution 
from AGN. The decrease is shown by a vertical line attached to a short horizontal line, and is 20%, on average, in both panels. The 
AGN luminosity used to correct the IR based SFR was determined using the relation derived in Murphy et al. (2009) between the AGN 
contribution to the IR luminosity and the difference between best-fit IR luminosities and those derived from 24 //m photometry alone (see 
§3.2.3). Sources that are detected in the hard band (2.0 — 8.0 keV) X-rays are identified by crosses. 
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Fig. 15. — The ratio of IR derived SFRs estimated from SED fitting (plus UV derived SFRs uncorrected for extinction) to extinction 
corrected UV SFRs against the color excess E{B-V) inferred from fitting the rest-frame UV (1250 - 2600 A) slope using UBViz optical 
imaging data for only those spectroscopically detected galaxies lying within a redshift range between 1.4 < z < 2.6 (i.e., approximately the 
BzK selection range). This is done to illustrate the difference in extinctions derived by using color relations (left panel; i.e., Daddi ct al. 
(2004)) to that when a proper slope is fit (right panel). The plus symbols in the right panel were generated using our revised extinction 
corrections (see §4). 



SFR by a factor of ^1.5. Comparison of the two methods as a function of redshift shows that the dispersion between 
them increases at the extremes of the redshift range for BzK selection (i.e., near z — 1.4 and 2.6). For comparison, 
the sources in the right panel of Figure 15 are re-plotted using the revised extinction corrections (i.e.. Equation 5) 
as plus symbols. Accordingly, using this color relation, rather than directly fitting the UV continuum for computing 
extinction corrections in this redshift range, does appear to introduce a systematic uncertainty for computing the UV 
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Fig. 16. — A compajrison of extinction estimates for galaxies lying within a redshift range between 1.4 < z < 2.6 (i.e., approximately the 
BzK selection range). Extinctions calculated using color relations from Daddi et al. (2004) are plotted against estimates from directly 
fitting the UV slope (see Figure 13). A one-to-one line is plotted as a dotted line. The direct fitting methods tends to yield larger reddening 
values. 



corrected SFRs. 
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